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There has been renewed theoretical interest recently in the Rabi model due to Braak's analytical 
solution and introduction of a new criterion for integrability. We focus not on the integrability of 
the system but rather why it is solvable in the first place. We show that the Rabi model is the 
limiting case of a more general finite dimensional system by use of a contraction and suggest that 
this is the reason for it's solvability, which still applies in the case of non-integrable but solvable 
variations. 



I. INTRODUCTION 

In classical mechanics the terms intcgrablc and solvable 
are often used interchangably, implying that a system 
that is integrable is also said to be solvable while a non- 
integrable system exibits chaos. In this sense it became 
clear in the early twentieth century that most classical 
dynamical systems are not analytically solvable, for ex- 
ample the perturbation series for the three body problem 
is only convergent in certain regions of the phase space. 
Later work by Kolmogorov, Arnold and Moser [1] estab- 
lished that this region covers a large volume for small per- 
turbations but has a complicated fractal structure. Thus, 
chaos and instability are still possible for small bodies in 
nearly Keplerian orbits and the solar system appears sta- 
ble because such bodies were either kicked out or fell into 
the Sun or Jupiter. 

An analogous understanding still does not exist in the 
case of quantum mechanics. The superficial observation 
that the problem to be solved is linear (the Schrodinger 
equation) misses the point that the Hilbcrt spaces of most 
systems of interest are infinite dimensional: a linear prob- 
lem in infinite dimensions has many of the analytical sub- 
tleties of finite dimensional non-linear ordinary differen- 
tial equations [2]. One way to approach the question of 
solvability of a quantum system is to ask if it has sensible 
finite dimensional truncations which are exactly solvable. 

Recently Braak [3] analytically solved the Rabi model 
[4], with Hamiltonian 



H, 



u>a'a 



-cr 3 + gai (a 



(1) 



which can be thought of as a two level atom with energy 
levels ±A/2 coupled to a quantized field with frequency 
lo. As Braak mentions, the general opinion was that this 
system is not solvable because the Hamiltonian is its 
only conserved quantity. A second conserved quantity 
would allow the Hamiltonian to be brought to diagonal 
form. For example, [5] in the Jayne-Cummings model 
Q = a^a + i<73 is conserved. In the Rabi model, al- 
though Q is not conserved, it only changes by a bounded 



amount: |AQ| < 1. We can see this by defining | q) as a 
basis in which Q is block diagonal, so that (q'\H\q) — 
if AQ = \q — q'\ > 1. This makes the Hamiltonian block- 
tridiagonal and we solve it in terms of continued frac- 
tions. Unlike Braak, who was interested in the integra- 
bility of the system, we ask another fundamental question 
"What makes a quantum mechanical system solvable?" 

To view this problem from a different perspective we 
look at the Hamiltonian 



H L = ujL 3 +AR s +gL 1 R 1 , 



(2) 



where L and R are angular momentum matrices with 
magnitude I and r respectively and we have absorbed a 
factor of 2 into g. This looks eerily similar to the Rabi 
model and in section V we use a contraction of the al- 
gebra to show that in the limit as I — > oo with r = -i, 
(2) becomes the Rabi Hamiltonian (1). The advantage 
of studying this new, more general Hamiltonian, is that 
unlike the Rabi model it has finite dimension. 

To find the matrix elements of Hl we use the basis 
states 



ifj L ) =| 1,7111, r,m r ) =\ l,mi)c 
-I, —I + 1, I and m r 



r, m ri 



(3) 



where mi = —I, —I + 1, I and m r = — r, — r + 1, r are 
the azimuthal components of angular momenta in the L3, 
R3 basis. This allows us to write the Hamiltonian as a 
block tridiagonal matrix 
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with 



A k = kuil r + AR 3 



B k = y/ 1(1 + 1) -k(k- TjgRi, 
where I r is an r dimensional identity matrix. 



II. CONTINUED FRACTIONS AND 
TRIDIAGONAL MATRICES 
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The general theory of analytical continued fractions 
was developed by Sticltjes in the late 19 th century while 
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studying divergent power series [6]. These continued frac- 
tions are of the form 



/(*) = 



1 



Z + di 



6? 



(6) 



Z + CL2 — - 



and are intimately connected with ordinary tridiagonal 
matrices of the form 



A = 



/ a bi 

ci oi b 2 

c 2 a 2 b 3 

c 3 a 3 

v . . . . 
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This connection can be seen by defining 

A_i(z) = l 
A (z) =a - z 

A fc (z) = (a k - z) A fe _i(z) - b k c k A k - 2 (z) 
which gives us a continued fraction 
A fe (z) 



(8) 



5 fc (z) 



A fc _iO) 
(a* - z) - 



b k c k 
Sk-i(z) 



whose roots are the eigenvalues of the matrix. If how- 
ever, A is an infinite dimensional matrix, the roots of 
S k represent the k th approximation to the eigenvalues of 
A. When a k are real and b k c k > the function S k is a 
Sturm sequence, meaning the zeros are real and the roots 
of S k -i arc the poles of Sk- This gives us an easy way to 
calculate higher approximations to the eigenvalues: be- 
tween any two pairs of poles of Sk is an eigenvalue. 

This tells us that as long as Sk is a convergent contin- 
ued fraction, even if the matrix it represents is infinite 
dimensional, it can still be solved to the desired level 
of precision. More importantly since the convergence 
of such continued fractions has been well established for 
over a hundred years this method gives a true check as to 
whether or not diagonalizing increasingly larger matrices 
will converge to the eigenvalues of an infinite dimensional 
matrix [6]. 



III. EIGENVALUES OF BLOCK TRIDIAGONAL 
MATRICES 

The usefulness of tridiagonal form to prove conver- 
gence and the block tridiagonal form of Hl prompts us 
to ask "is it also possible to develop relations similar 
to (8) and (9) for block tridiagonal matrices?" The an- 
swer is, in some cases, yes. Using the transfer matrix 
method of Molinari [7] we can find the eigenvalues of an 



(n + 1) X (n + 1) block tridiagonal matrix 

\ 
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where A k ,B k are mx m square matrices, and A k is of 
the form A k — zl m , where z are the eigenvalues. In this 
case detM = 0, allowing us to set 



/ A Bx 

' Ci At B 2 

C 2 A 2 B 3 

C 3 A 3 

V . . . . 
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(11) 



where f is a null vector with components ip k G C", 
giving us a set of equations 

Aoij + Bi^i = 
B k+1 4' k+ i + A k 4' k + Cktpk-i = 

A n ip n + C n il) n+1 = 0. (12) 

This can be written recursively as 



(9) 


ipk+i 






ipk 





-B k+1 A k -B^Ck 





which defines the transfer matrix 



ipk-i 



(13) 



T k = 
T n = 



1 1)1 



Ira 







—B n A n _x —B n x Cn- 

Ira 



-B^A -B^ 

Ira. 



If we define the top left element of T k as 

Tk,u = —B k+1 AkTk-i,ii — B^ +1 CkTk-2,u, 
since ip n +2 = an d ip-i = 





lpn+1 



Tpl 





(14) 



(15) 



(16) 



so detT„ ; ii = is the same as detM = 0, defining the 
eigenvalue equation. 



IV. SPECTRUM OF H L 

For the most general form of Hl, integer values of r 
correspond to singular matrices that have no inverse, 
so the above method can be modified to the form of 
Salkuyeh[8]. It is possible to solve Hl for any half integer 
spin but the most elegant case, which is also interesting 
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because of it's ties to the Rabi model Hr, is for r = |. 
In this case 



Ri 
A, 



A 

kl + -C73 



C k = 5a /Z(Z + 1) - - lK 



D 



(17) 

where we have absorbed a factor of ^ into g and defined 

bk = g + 1) — k(k — 1). This gives us the eigenvalue 
equation dct (Hl — zl) = = det (Tu(z)) , where the 
transfer matrix can be simplified by multiplying each 2x2 
matrix in Tk by bk to give 



T, 



Tk-i 



oi (A k - zl) b k I 
-b k 

Tk,ii = (T\AkTk-\,\\ — b 2 JTk-2,\\ 



(kuj - z) crj - iy(T 2 ) 7fe-i,ii 



Due to the special property of Pauli matrices OiOj 
Sij + ieijk&k we see that if we define 

S-m — Jm , 1 1 — 1 1 1 

we have the recursive matrix 



S m = (mcj — z)(7i — iAcr 2 — BfS, : 1 



m-li 



a(z) 
b{z) 



which is of the form 
with inverse 

This allows us to write S m in matrix form 
/ a+ 




(18) 
T J = 

(19) 
(20) 

(21) 
(22) 



where 







(23) 

muj — z ± — , and by setting S m = we 



obtain two terminating continued fractions 

A bl 

2 S, 



(z) = muj — 2=F 



m— 1,± 



(24) 



where m = + whose zeros are the eigen- 

values of the Hl. Because these are finite dimensional 
matrices, calculating the roots of Si _|_ gives the even par- 
ity spectrum while the roots of S^— give the odd parity 
spectrum. 



V. CONTRACTION OF H L TO H R 

We could of course perform the same procedure on the 
Rabi model Hr to find it's eigenvalues in a similar way 
but a deeper connection between the two systems can be 
seen by performing a singular change of basis. It was 
discovered by Inonu and Wigner [9] that a transforma- 
tion of this type changes one Lie Algebra into another 
using a process called contraction. We briefly summarize 
Gilmorc's[10] description of a different class of contrac- 
tions to show the relationship between (1) and (2). 

A Lie algebra defined by the basis vectors Xi is closed 
under commutation, so the commutators 



[Xi,Xj] - C^Xk 



(25) 



are contained in the algebra. The structure constants 
completely determine the algebra, however it is possible 
to perform a change of basis transformation 



Y, = M/Xi 



where the new structure constant C,-,- becomes 



a 



(M 



(26) 



(27) 



due to a non-singular transformation. If we allow the 
transformation to be parameter dependent, where 



Y i =M?{e)X j 



C 



(28) 



the structure constant often converges to a new Lie Al- 
gebra if (e) becomes singular in the limit as e — > oo. 

One representation of angular momentum is the com- 
pact unitary group U(2), which is spanned by the oper- 
ators J 3 , J±, J , (Jo is the identity) which correspond to 
the commutation relations 

[J 3 , J±] = ±J± 
[J+, J_] = 2J 3 
[Jo,J] = 0. (29) 

We can now change basis to the Heisenberg group H4 by 
using 



(30) 









(£). 




( J l) 


\ho) 




J 3 
[Jo/ 



] = ±h± 

2c 2 h?. 



which gives us 

[h 3 , h± 
[h + ,h j 
[h ,h]=0 

and in the limit as c — > 

[ha,h±] = ±h± 
[h+,h-\ = -h - 



(31) 



(32) 
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which satisfy the same commutations as the single mode 
photon operators 

\N = at a, a] = —a 
[N = a! a, at] = a* 

[a^a] = -1. (33) 

We can now identify (in the limit as c — > 0) 



So a contraction on Hl changes (24) to 



h 3 = N 
h + = at 
h- = a. 



(34) 



To see how the basis states change we first operate h% on 
the angular momentum state | j, m) to get 

^3 I j, m) = y J 3 + ^ J ^ I j, m) = + | j, m). 

(35) 

The ground state of this system corresponds to m = —j , 
so the n th state is n = j + in. In order for the limit to 
be well defined we require 



lim I m + — !- j = Hm ( n — j + — ^ 

c^o V 2c 2 / c^o V 2c 2 



(36) 



to also be well defined. We have already equated /13 with 
the number operator, so the requirement becomes 



(37) 



telling us that 2jc 2 = 1, or in other words as c — > 0, 
j — > 00, and 



^,±(2) = few - z =F (-1 



2 5 fc „!,±(z) : 



(40) 



fc = 0,1,2,3,... Its zeros give us the spectrum of the 
Rabi model Hr. As discussed in section II, there is 
a rapidly convergent algorithm to find the zeros, which 
takes advantage of the fact that there is a zero of Sk,±(z) 
in between two of its poles, which are simply the zeroes 
of the previous approximation, Sk-i,±(z). Thus we can 
limit the search for each zero to these intervals, increasing 
the size k of the matrix in each step. 

VI. CONCLUSION 



We show that the Rabi model can be thought of as the 
limit of a sequence of finite dimensional block tridiagonal 
Hamiltonians, each of which can be solved by a contin- 
ued fraction method. The solvability of the Rabi model 
can thus be understood as due to this finite dimensional 
truncation and the existence of approximate conserva- 
tion laws (selection rules for transition matrix elements) 
that ensure block tridiagonality, explaining why even a 
Hamiltonian with a broken symmetry can be solved [3]. 
Conversely, we should expect that systems which do not 
allow convergent finite dimensional approximations ex- 
hibit quantum chaos. We hope to construct such an ex- 
ample in a later publication. 



lim h 3 \j,m)=n\ 00, n). 
c— >o 



(38) 



Similarly, we see that 

at I n) — lim cj + I j, m) 

c— >o 



lim \f j{j + 1) — m(m + 1) | j, m + 1) 

c— >o 



lim J(l - c 2 n)(n + 1) I j, m + 1) 



Vn + 1 I n + 1). 



(39) 
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